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ABSTRACT 

We analyze spectra for a large sample of 578 Active Galactic Nuclei to examine the relationships 
between broad emission line properties and central supermassive black hole (SMBH) mass. We estimate 
SMBH masses by applying the virial theorem to the C IV A1549 broad emission line. Although the 
FWHMs of C IV and HP appear nearly unrelated in individual objects, these FWHMs are well correlated 
when averaged over sub-samples in our database. Therefore, the lines are equally valid indicators of the 
average SMBH mass in quasar samples. Our sample spans five orders of magnitude in SMBH mass, 
six orders of magnitude in luminosity, and a redshift range from < z < 5. Most lines diminish in 
equivalent width with increasing black hole mass (the usual "Baldwin Effect" ) and there are no trends 
with redshift. Recent studies indicate that there is a relationship between SMBH mass and the overall 
bulge/spheroidal component mass of the surrounding galaxy. This relation, together with the well-known 
mass-metallicity relationship among galaxies, predicts a relationship between SMBH mass and quasar 
metallicity. We estimate the metallicity in the broad emission line region by comparing several line ratios 
involving nitrogen to theoretical predictions. We find that the data are consistent with a trend between 
SMBH mass and metallicity, with some line ratios indicating a very strong trend, but the uncertainties 
in several other important line ratios are too large to confirm or test this correlation. 

Subject headings: galaxies: active — quasars: emission lines — galaxies: formation 



1. INTRODUCTION 

Supermassive black holes (SMBHs) are believed to drive 
the tremendous energy output of quasars, and more gen- 
erally, active galactic nuclei (AGNs). Recent studies have 
shown that the centers of all massive galaxies today har- 
bor SMBHs (Gebhardt et al. 2000a; Ferrarese & Merritt 
2002). Moreover, the mass of the black hole scales di- 
rectly with the mass of the host galaxy, or more specif- 
ically, the mass of the host galaxy's bulge or spheroidal 
component. The early results were based on a correlation 
between SMBH mass and bulge luminosity, Lbuigc (Magor- 
rian et al. 1998; Laor 1998; Wandel 1999), which showed a 
large scatter, as much as two orders of magnitude in SMBH 
mass (Ferrarese & Merritt 2000). Subsequent studies fo- 
cused on the relationship between the masses of SMBHs 
and the velocity dispersions, ti, of their host bulges or 
spheroids. This relation shows considerably less scatter 
than the Msmbh — ^buigc relation (Ferrarese & Merritt 
2000; Gebhardt et al. 2000a; Merritt & Ferrarese 2001; 
Tremaine et al. 2002). There is also no significant dif- 
ference between the relationships measured for active and 
quiescent galaxies (Gebhardt et al. 2000b; Ferrarese et 
al. 2001). Recently, SMBH mass has also been shown to 
correlate strongly with the global structure of ellipticals 
and bulges (e.g. more centrally concentrated bulges have 
more massive SMBHs). This relationship is as strong as 
the Msmbh — cr relationship with comparable scatter (Gra- 
ham et al. 2001; Erwin et al. 2002). Studies that carefully 
model the bulge light profiles of disk galaxies and thereby 
obtain more accurate values of ibuige also show less scat- 
ter in Msmbh — ^buigc, similar to that in the Msmbh — o- 
relationship (McLure & Dunlop 2002; Erwin et al. 2002; 



Bettoni et al. 2003). 

The overall mass of a galaxy is also known to corre- 
late with its metallicity (Faber 1973; Zaritsky et al. 1994; 
Jablonka et al. 1996; Trager et al. 2000). This relation- 
ship is generally understood in terms of the gravitational 
binding energy. More massive galaxies have deeper po- 
tential wells, making it harder for supernovae and stel- 
lar winds to blow matter out of the region. Therefore, 
the gas remains in the galaxy longer, where it is fur- 
ther reprocessed by stars, leading to higher metallicities. 
This mass-metallicity relationship, combined with the re- 
lationship between SMBH mass and galaxy mass, suggests 
that there should be a relationship between SMBH mass 
and the metallicity of the gas clouds surrounding quasars 
(Hamann & Ferland 1993). 

The metallicity of the gas in the broad emission line 
region (BLR) of quasars can be estimated by analyzing 
prominent emission lines. BLR metallicities should be rep- 
resentative of the gas in the central regions of galaxies 
if the BLR was enriched by stars in those environments. 
Previous studies have demonstrated that we cannot sim- 
ply examine ratios of strong metal to hydrogen lines, such 
as C IV A1549 / Lya, to measure the C/H abundance. 
This is because strong collisionally excited lines like C IV 
play an important role in cooling the gas. Increasing the 
metallicity and thus C/H leads to lower gas temperatures 
and an almost constant C IV / Lya ratio (Hamann & 
Ferland 1999). The best abundance diagnostics for the 
BLR involve ratios of lines that are i) not important in 
the cooling, and ii) have similar excitation/emission re- 
quirements (Hamann et al. 2002). Ratios of emission lines 
involving nitrogen, N, are especially valuable in determin- 
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ing the metallicity, Z, because of the expected "secondary" 
N production via the CNO cycle of nucleosynthesis in stars 
(Shields 1976; Hamann & Ferland 1992, 1993, 1999; Fer- 
land et al. 1996; Hamann et al. 2002). In the CNO cy- 
cle, nitrogen is produced from existing carbon and oxygen. 
The nitrogen abundance should therefore scale roughly as 
N/H oc Z2 or N/0 oc 0/H oc Z (Tinsley 1980), providing a 
sensitive metallicity diagnostic even when direct measures 
of Z = O /H are not available. Observations of H II regions 
indicate that secondary nitrogen production and N/O oc 
0/H scaling dominate for Z > 0.2 — 0.3 (van Zee et al. 
1998; Pettini et al. 2002). 

Some other studies suggest that there can be significant 
departures from the simple N/O oc 0/H relationship ow- 
ing to time-dependent or perhaps metallicity-dependent 
yields of N, O, and C (eg. Henry et al. 2000). These 
effects are often invoked to explain the observed object- 
to-object scatter in N/O versus 0/H. However, there is a 
general consensus that large (solar or higher) N/O ratios 
indicate large (solar or higher) metallicitics. Moreover, the 
most homogeneous and therefore most reliable H II region 
data (Pilyugin et al. 2002; Pettini et al. 2002) indicate 
that the scatter in N/O decreases and the validity of N/O 
oc 0/H improves in the regime of roughly solar or higher 
metallicities. This is the regime of the gas near quasars, 
based on studies of both their broad emission and intrinsic 
narrow absorption lines (Constantin et al. 2001; Dietrich 
ct al. 1999, 2003a; Hamann 1997; Hamann & Ferland 
1999; Hamann et al. 2002; Osmer et al. 1994; Petitjean 
et al. 1994; Warner et al. 2002). In this paper, we simply 
assume that N/O oc 0/H applies. We also note that our 
analysis of composite quasar spectra naturally averages 
over object-to-object variations. Any metallicity trends 
that we infer from the emission line ratios depend on this 
assumption of N/O oc 0/H in addition to measurement 
and photoionization modeling uncertainties. 

We have collected spectra of over 800 "Type 1" AGNs 
(quasars and Seyfert 1 galaxies) from the HST archives 
and various ground-based observing programs to exam- 
ine trends in their emission line properties (Dietrich et al. 
2002). 578 of these spectra span the rest-frame UV wave- 
lengths needed for this study. We compute five composite 
quasar spectra from this sample, representing different in- 
tervals in SMBH mass from 10*^ to 10^° M©. We present 
measurements of the emission lines in these composite 
spectra and investigate the relationship between SMBH 
mass and the metallicity of the BLR gas. 

2. SMBH MASS DETERMINATIONS 

We estimated black hole masses by applying the virial 
theorem, Msmbh — rv^ jG, to the hue-emitting gas, where 
V is the velocity dispersion and r is the radial distance away 
from the SMBH (Peterson & Wandel 2000). The relation- 
ship between v and the FWHM of the broad emission line 
profile depends on the kinematics and geometry, but can 
be expressed as v = fc x FWHM, where fc is a factor on 
the order of unity (McLure & Dimlop 2001: Vcstcrgaard 
2002). Kaspi et al. (2000) express the SMBH mass as 



M = 1.464 X lO-^M, 
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FWHM applies to the broad emission line profile, and 
^?BLR is the radial distance between the BLR and the cen- 



tral mass/continuum source. Equation (1) assumes the gas 
in the BLR is gravitationally bound and BLR velocities are 
random. The exact geometry and kinematics are still de- 
bated (e.g. McLure & Dunlop 2001), but are unknown, so 
we take this approach of random velocities (fc = ■v/3/2) be- 
cause it is the most basic and allows for comparison with 
most other work. Any offset in mass due to different geo- 
metric or kinematic assumptions would only affect the ab- 
solute SMBH masses and not any trends derived between 
metallicity and SMBH mass. We estimate i?BLR based on 
the observed relation between i?BLR for a particular line 
and the continuum luminosity (Kaspi et al. 2000; Vester- 
gaard 2002; Wandcl ct al. 1999). A particular line must 
be specified because reverberation studies have shown that 
the BLR is radially stratified, such that higher ionization 
lines tend to form closer to the central engine than lower 
ionization lines (Peterson 1993). Kaspi et al. (2000) find 
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where Rbv,k{H(3) is the radial distance between the cen- 
tral continuum source and the H(3 emission line region 
and AiA(5100A) is the continuum luminosity at 5100 A 
in the AGN rest frame. We choose to use the C IV A1549 
emission line instead of H(5 to estimate the SMBH mass 
because it is more readily observed across the entire red- 
shift range from z~0toz~5. CIV may in fact yield a 
better estimate of SMBH mass than Hf3 because H(3 can 
include a narrow line component that forms in a different 
region and thus must be subtracted before measuring the 
FWHM. C IV has no such narrow component (Wills et al. 
1993, Vestergaard 2002). It has recently been suggested 
that Mg II A2798 may be a better indicator of SMBH mass 
than C IV (McLure & Jarvis 2002). However, C IV is ob- 
servable over a wider range of redshifts and is clearly the 
best choice among UV lines at shorter wavelengths. Sig- 
nificant blending effects are much less common for C IV 
than other strong lines in the rest-frame UV spectrum and 
absorption effects are much less common for C IV than for 
Lya (Vestergaard 2002). 

Reverberation studies indicate that the size of the broad 
emission line region for C IV is about half that of H[3 
(Stirpe et al. 1994; Korista et al. 1995; Peterson 1997; 
Peterson & Wandel 1999). We modify Equation (2) to ac- 
count for this and we use a powerlaw of the form Fn oc j^" 
with a = —0.4 to estimate a conversion between the 
continuum luminosity at 5100 A and 1450 A. We select 
a = —0.4 based on average quasar spectra from Brother- 
ton et al. (2001), Vanden Berk et al. (2001), and Dietrich 
et al. (2002). This yields 

ALa(1450A) 



i?BLR(CIV) = 9.7 



n 0.7 



It - days (3) 



.10^4 ergs s- 

Vestergaard (2002) finds that this technique of using the 
continuum luminosity and the FWHM of C IV yields an 
estimate of the central black hole mass that has a Icr uncer- 
tainty of a factor of three when compared to reverberation 
mapping studies that use Hfi to estimate SMBH masses. 
From Equations (1) and (3), we obtain 

/ FWHM(CIV) Y , 

V 103 km s-i J 1^1044 ergs s"! ) 
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This is essentially the same as (within 10% of) the 
mass relationship derived by Vestergaard using C IV and 
ALa(1350A) (For this comparison, we again assTime the 
powerlaw index, a. = —0.4 between 1450 and 1350 A.) The 
difference of ^ 10% between the two formiilae is negligible 
compared to the la uncertainty of a factor of three. See 
also Netzer (2003) and Corbett et al. (2003) for further 
discussion of the uncertainties. 

While multi-epoch spectra are preferable (to average 
over variabilities) , single-epoch spectra may be used if the 
signal-to-noise ratio is high enough because the FWHM 
of single-epoch spectra are generally consistent with the 
FWHMs of the mean and rms spectra obtained from multi- 
epoch observations (Vestergaard 2002). It is important 
to accurately measure the FWHM because the uncertain- 
ties in SMBH mass are dominated by the uncertainties in 
the FWHM. Our composite spectra are each made up of 
single-epoch observations of multiple quasars and there- 
fore should be approximately equivalent to multi-epoch 
observations of the same sources. 

3. DATA & ANALYSIS 

We have compiled a sample of more than 800 quasar 
spectra, spanning a redshift range from < z < 5 and 
six orders of magnitude in intrinsic luminosity. The spec- 
tra were obtained by several groups using various ground- 
based instruments, as well as the International Ultraviolet 
Explorer (lUE) and the Hubble Space Telescope (HST). 
One unique aspect of this sample is that it contains new 
observations of faint quasars at redshift z > 2.5 (Diet- 
rich et al. 2002). Thus we can avoid to some degree the 
bias toward higher luminosities at higher redshifts that af- 
fects other magnitude-limited samples. 578 of the spectra 
have UV wavelength coverage that encompasses the range 
950 < A < 2050 A. Each spectrum was transformed to the 
rest-frame using a redshift measured from the centroid of 
the upper 50% of the C IV A1549 profile. See Dietrich et 
al. (2002) and Dietrich et al. (2003, in prep) for more 
details on the AGN sample and how luminosities were ob- 
tained. Throughout this paper, we use the cosmological 
parameters Hq = 65 km s~^ Mpc~^, Om = 0.3, and Oa = 
(Carroll, Press, & Turner 1992). We used radio flux den- 
sities given in Veron-Cetty & Veron (2001) to determine 
the radio loudness for the quasars. 

We use an automated program to estimate the FWHM 
of C IV in each spectrum. The program averages over 
noise by applying a smoothing routine that calculates the 
median of each five-pixel interval. It then fits a local lin- 
ear continuum (contrained by the flux in 25 A wide in- 
tervals centered around 1442.5 A and 1712.5 A) around 
C IV, measures the peak flux of the emission line (of 
the smoothed spectrum), and uses those to estimate the 
FWHM of the fine. Comparisons between the FWHMs es- 
timated by the program and those measured manually in- 
dicate an error of < 10% in the estimates. Lines containing 
significant absorption are flagged by the program and their 
FWHMs are estimated manually (by manually interpolat- 
ing across the absorption feature). We use the FWHM of 
C IV and the continuum luminosity, ALa(1450A) to es- 
timate the central SMBH mass of each quasar based on 

^ The Image Reduction and Analysis Facility (IRAF) is distributed 
the Association of Universities for Research in Astronomy, Inc. (USI 



Equation (4). 

We then sorted the quasars by SMBH mass into five 
bins: 10^-10^ M0, IO^-IO^Mq, IO^-IO^Mq, 
Mq, and > 10^° M©, and computed five composite spec- 
tra. Each composite spectrum is the average of all the 
quasar spectra in a bin (see Table 1 for data about each 
composite spectrum, including the powerlaw index, a, of 
the continuum fit, the number of objects contributing to 
the composite at and FWHM of C IV and i?/3, continuum 
luminosity, and mean SMBH mass). Calculating compos- 
ite spectra increases the signal-to-noise ratio significantly, 
making it easier to measure weaker emission lines that may 
be lost in the noise in individual spectra. 

Since narrow absorption features may influence the 
emission line profiles in composite spectra, we developed a 
method to detect strong narrow absorption features that 
exclude the contaminated spectral region of individual 
spectra from the calculation of the composite spectrum. 
The procedure is roughly as follows (see Dietrich et. al 
2002 for more details). First, a preliminary mean spec- 
trum is calculated from all of the spectra in a given mass 
bin. Each individual spectrum is then divided by this pre- 
liminary mean to derive a ratio spectrum. The ratio spec- 
trum is smoothed with a running boxcar function, which 
provides in addition to the smoothed ratio also a root- 
mean-square (rms) value for each wavelength box. A nar- 
row absorption line will cause a sharp increase in the rms 
spectrum. The contaminated spectral regions of individ- 
ual spectra, identified by these rms spikes, are excluded 
from the calculation of the final composite spectra, and 
no interpolation has been applied (Dietrich et al. 2002). 

We next corrected each composite spectrum for strong 
iron emission lines. To get a first estimate of the contribu- 
tion of Fe emission, we used the empirical Fe emission tem- 
plate that was carefully extracted from IZwl by Vester- 
gaard & Wilkes (2001), which they very kindly provided 
for this study. The template is mostly Fe II omission but 
contains some Fe III emission as well. Wc have modified 
this Fe emission template int he UV range AA1200— 1540 A 
in accordance with photoionization model calculations pre- 
sented by Verner et al. (1999). Based on our analysis of 
quasar composite spectra (Dietrich et al. 2002), we found 
that the strength of Fe line emission in the wavelength 
range AA1420 - 1470 A is equal to ~2 % of the integrated 
pseudo continuum flux in this wavelength range (Dietrich 
et al. 2002). Fe contributions vary for other wavelength 
ranges, but we scale the entire template to match the 2% 
contribution in the range 1420 - 1470 A. This estimate of 
a 2 % contribution is conservatively low to avoid over- 
subtraction of the Fe emission template. The spectral 
width of the Fe emission features was adjusted to a width 
corresponding to the FWHM of the C IV A 1549 emission 
line profile of each quasar. The appropriated scaled Fe 
emission template was subtracted from the corresponding 
composite spectrum. In spite of its small contribution, 
the estimate of the Fe line emission improves particularly 
the measurements of N III] A1750 and lines such as He II 
A1640 near an unidentified A1600 A emission feature (Laor 
et al. 1994; Vestergaard & Wilkes 2001). 

We used the task NFITID in the IRAF^ software pack- 

the National Optical Astronomy Observatories, which is operated by 
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age to fit the continuum of each Fe-subtracted spectrum 
with a powerlaw of the form oc v". The fits were con- 
strained by the flux in wavelength intervals between the 
emission lines, namely in 25 A wide windows centered at 
1460 A, 1770 A, and 2000 A. The spectral indices overall 
differ from a = —0.4 because the powerlaw fit is to a nar- 
rower wavelength range than, for example, Dietrich et al. 
(2002) or Vanden Berk et al. (2001). However, there is 
a significant trend for steeper (softer) UV spectra in the 
low-mass sources (see Table 1). Figures 1 and 2 show the 
final Fe Il-subtracted composite spectra normalized by the 
continuum fits. 

To measure the broad emission lines, we use the spec- 
tral fitting routine SPECFIT (Kriss 1994), which employs 
minimization. We fit each line with one or more Gaus- 
sian profiles, with the understanding that each individual 
Gaussian component may have no physical meaning by 
itself. Our goal in fitting the lines is simply to measure 
the total line strengths free of blends. Our strategy is 
to obtain the best fit with the least number of free para- 
maters and, when necessary, to use the profile of strong un- 
blended lines, such as C IV, to constrain the fits to weaker 
or more blended lines. Figure 3 shows an example of these 
fits. The Appendix provides details of our fitting proce- 
dure for different lines/blends. Table 2 lists the resulting 
line fluxes relative to Lya, plus the rest-frame equivalent 
widths (REWs) as measured above the fitted continuum, 
and the FWHM of each line. The fluxes, REWs, and 
FWHMs given in Table 2 are from the total fitted pro- 
files, which can include several multiplet components and 
up to three Gaussian profiles per component (as described 
in the Appendix). 

The primary uncertainty in our flux measurements is 
the continuum location. We estimate the 1 a standard de- 
viation of our measurements of the flirxes of C IV and Lya 
to be <10% based on repeated estimates with the contin- 
uum drawn at different levels. By the same method, we 
estimate the uncertainty in N V, N III], He II, O III], and 
other weaker lines to be ~10-20%. These estimates do not 
include the uncertainties due to line blending, which can 
be important for some of the weak lines and for N V in 
the wing of Lya. The blending and continuum placement 
uncertainties are both greater for the higher mass spectra 
because of the broader profiles and generally lower REWs. 

4. C IV VS. H/3 AS A MASS DIAGNOSTIC 

To test the hypothesis that i?BLR for C IV is typically 
half the value of that for H(}, we compare the measured 
FWHM of C IV to that of F/3 (see Figure 4). If the ve- 
locities are virialized and i?BLR for C IV is half that of 
Hj3, we expect the FWHM of C IV to be a factor of ^2 
larger than the FWHM of Hp. FWHMs are estimated by 
an automated program as described in §3. There can be 
separate BLR and Narrow Line Region (NLR) contribu- 
tions to (Marziani et al. 1996, Hamann ct al. 1997), so 
each Hp line is looked at individually. In 5 out of 74 cases, 
where there is obvious narrow HP emission with the same 
redshift and profile as [O III] A5007, this narrow emission 
is clipped off and the FWHM of HP is manually estimated. 
There appears to be no real correlation between the two 
FWHMs for individual objects, but averages of objects in 
selected ranges of C IV FWHM fall generally between a 1:1 



correlation and the expected \/2:l correlation (see Figure 
4). It should be noted that observations of HP and C IV 
for individual objects were not necessarily made simulta- 
neously, which may contribute to the scatter in Figure 4. 
Five composite spectra, created from objects in selected 
ranges in SMBH mass (derived from C IV using Eqn. 4 
and plotted as open boxes with error bars), do appear to be 
consistent with the expectation that FWHM(C IV) = \/2 
FWHM(i//3). All spectra contributing to the composite 
spectra contain C IV, but only some contain HP (see Table 
1). The narrow component of Hp was not subtracted from 
any spectrum before the composite spectra were created. 
In the composites, the narrow HP contribution smoothly 
blends into the broad-line profile. Filled diamonds rep- 
resent the averages of objects in selected ranges of C IV 
FWHM. These averages contain only the individual ob- 
jects (plotted with -I- signs in the figure) that have data 
at both C IV and Hp. Thus, while there is a large scatter 
among individual sources, the FWHMs of C IV and Hp 
generally scale as expected in both the composite spectra 
and the averages of individual objects. Also note that any 
offsets in Fig. 4 affect only the absolute SMBH masses, 
with an uncertainty of up to a factor of two, but do not 
affect any trends with SMBH mass. 

We also estimate the SMBH mass derived from C IV 
and the SMBH mass derived from Hp. Masses are de- 
rived from Hp using Eqns. 1 and 2. We compare the two 
masses for all objects in our sample that have observations 
of both lines (see Figure 5). Filled diamonds represent the 
averages of objects in selected ranges of SMBH mass (as 
derived from C IV). As in Fig. 4, these averages con- 
tain only the individual objects that have data at both 
C IV and Hp. We find that that there is approximately 
a 1:1 correlation between the mass obtained from C IV 
and that obtained from Hp. There can be significant de- 
viations from this for individual objects, but the relation 
holds well for averages of many objects. Hereafter, we will 
discuss only the masses obtained from C IV. 

5. RESULTS 

5.1. SMBH Mass vs. Redshift, Luminosity, and FWHM 

Figure 6 shows the redshift distribution of the entire 

sample as a function of SMBH mass. The sample was 
designed to span the widest possible range of luminosity 
at each redshift, and the distribution of points in Fig. 6 
closely follows the distribution in L vs. z (see Dietrich et 
al. 2002). There is no trend in SMBH mass with redshift. 
We expect SMBH mass to correlate with both continuum 
luminosity and FWHM because we derived our mass esti- 
mates from these quantities (Eqn. 4). Continuum lumi- 
nosity is plotted as a function of SMBH mass in Figure 
7, and shows the expected trend with little scatter. This 
relationship is much tighter than the relationship with the 
FWHM of C IV (see Figure 8), even though our SMBH 
mass estimates depend more sensitively on FWHM than 
luminosity. This may be due to the larger dynamic range 
in luminosity compared to FWHM (sec Netzer 2003 for 
further discussion). In all three plots (Figs. 6-8), no trend 
is seen based on radio loud vs. radio quiet objects. 

5.2. Line Properties versus SMBH Mass 
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The Baldwin Effect is an observed anti-correlation be- 
tween continuum luminosity (and thus SMBH mass, since 
that is estimated based on luminosity) and the equivalent 
widths of emission lines such as C IV, Lya, C III], and 
O VI (see Baldwin 1977; Groom et al. 2002; Dietrich et 
al. 2002; Espey et al. 1993; Osmer et al. 1994, Osmer 
& Shields 1999; Zheng & Malkan 1993). Interestingly, the 
Baldwin Effect is not seen in N V even though it seems to 
be present to some extent in N III]. Dietrich et al. (2002) 
disciiss the Baldwin Effect of this same sample sorted by 
continuum luminosity. The Baldwin Effect is also evident 
in our spectra sorted by Msmbh- Figure 9 plots rest-frame 
equivalent widths as a function of SMBH mass (see also 
Figs. 1 & 2). Al III curiously does not show the Bald- 
win Effect in these composite spectra and, if anything, 
its equivalent width actually increases with SMBH mass. 
This is contrary to its behavior in Dietrich et al. (2002) 
when sorted by continuum luminosity. The rest of the lines 
plotted in Fig. 9 show generally the same behavior as in 
Dietrich et al. (2002). Dietrich et al. (2002) find that 
the slope of the Baldwin Effect becomes steeper in most 
lines at higher luminosities. We observe a similar trend 
for steeper slopes at higher SMBH masses with the no- 
table exception of N III] , which slightly increases in REW 
at higher SMBH masses. 

5.3. Metallicity versus SMBH Mass 

Figure 10 shows the metallicities inferred by comparing 
emission line flux ratios for each composite spectrum to 
the theoretical results in Hamann et al. (2002). We ob- 
tain our metallicity estimates by comparing the emission 
line flux ratios to plots of metallicity vs. line ratio based 
on theoretical models (see Figure 5 in Hamann et. al 2002 
and Figure 3 in Warner et al. 2002). We prefer to estimate 
the metallicities based on the calculations in Hamann et 
al. that use a segmented powerlaw for the photoionizing 
continuum shape. This continuum shape is a good ap- 
proximation to the average observed continuum in quasars 
(Zheng et al. 1997; Laor et al. 1997). It also yields inter- 
mediate results for line ratios, such as N V/He II, that are 
sensitive to the continuum shape (see below). 

The uncertainties in the metallicities shown in Figure 
10 derive simply from the l-cr measurement uncertainties 
discussed in §3. They do not reflect the theoretical uncer- 
tainties in the technique we use to derive metallicities from 
the line ratios (see Hamann et al. 2002 for discussion). 
Also note that there is a larger statistical uncertainty in 
the lowest mass composite spectrum, because it is com- 
posed of only nine objects (the small number of spectra 
in this composite make it easier for one or two objects to 
dominate the composite). 

Our best estimate of the overall metallicity from each 
spectrum (labeled as Average in Fig. 10) is obtained by av- 
eraging the results of N III]/C III], N III]/0 III], N V/C IV, 
and N V/0 VI. We select these ratios because we believe 
they are the most accurately measured ratios and the most 
reliable from a theoretical viewpoint (Hamann et al. 2002). 

We believe the ratios involving N IV] arc less reliable 
for several reasons. First, the N IV] flux is difficult to 
measure because this line is weak and in the wing of the 
much stronger C IV line profile. Second, ratios such as 
N IV] /C IV, which compare a permitted and intercombi- 



nation line, are potentially more sensitive to uncertainties 
in the line optical depths and radiative transfer. Third, 
perhaps as a consequence of the greater uncertainties in 
N IV], there is more object-to-object scatter in the re- 
sults derived from N IV] compared to N III] (Hamann et 
al. 2002). We also exclude N V/He II from our estimate 
of the average metallicity because i) this ratio compares 
a coUisionally excited line (N V) to a recombination line 
(He II), and ii) it can yield (for this continuum shape) very 
high metallicities that are obtained by extrapolations be- 
yond the theoretical results calculated by Hamann et al. 
(2002). 

All of the line ratios involving N III] and N V show N/0 
and N/C ratios that are solar or greater. This implies a 
metallicity of > 1 Z© if N is mostly secondary. 

5.4. Do NIII] and N V Disagree? 

All four ratios involving N V show a strong trend in 
metallicity with SMBH mass. The slope of this trend is 
not linear and seems to increase as mass increases. The 
ratios involving N III], however, show no apparent correla- 
tion between metallicity and SMBH mass. We believe that 
this apparent discrepancy between the two ratios may be 
due to the large FWHMs and small REWs of the emission 
lines in the composite spectra representing higher SMBH 
masses. 

The spectra representing higher SMBH masses have 
broader lines because the mass derivation depends on the 
line FWHM (Eqn. 4). This may lead to systematic over- 
estimating of the continuum if the wings of the lines blend 
together and never quite reach the continuum level. In 
particular, notice the interval around 1460 A (see Figs. 1 
and 2) that we used to constrain our continuum fits (see 
§3). In the higher mass spectra, C IV and Si IV -|- O IV] 
seem to blend together, never quite reaching the true con- 
tinuum, and instead forming a "U" shape. This effect 
of overestimating the continuum could lead to underesti- 
mates of N III] and N IV] in the higher mass spectra and 
thus negate or weaken any real trend between SMBH mass 
and the flux ratios involving these lines. For example, we 
flnd that a drop in the continuum level of only 5% will 
more than double the measured flux in N III] and N IV], 
but increase the flux in stronger lines by only ~ 15% to 
30%. 

In order to more thoroughly examine the effects of 
FWHM on the measured line fluxes, we sort the quasar 
spectra by the FWHM of C IV and create composite spec- 
tra for different ranges in FWHM. We then fit the con- 
tinuum and measure the emission line fluxes as described 
in §3. The resulting line ratios are plotted in Figure 11. 
We find that the flux ratio of C III]/C IV remains roughly 
constant at different line widths, but N III]/N V declines 
as the lines get broader. Because of this, we again see a 
discrepancy in metallicity trends between ratios involving 
N III] and those involving N V. The ratios involving N V 
indicate increasing metallicity with increasing line width. 
This is expected because objects with broader emission 
lines should have higher black hole masses and thus higher 
metallicities (§1). The N III] ratios, however, show a "U"- 
like behavior in Figure 11, with the narrowest lines having 
the highest metallicity. The composite spectrum with the 
narrowest lines is also the only one in which the metal- 
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licities derived from N III] ratios agree with those derived 
from N V ratios. 

There are two possible explanations for the discrepancy 
between the line ratios. The first possibility is that N V 
is somehow enhanced in quasars with broader emission 
lines (and/or larger L and Msmbh) by a process that we 
do not understand. The second possibility is that we are 
systematically overestimating the continuum and thus un- 
derestimating the flux of N III] in spectra with broader 
lines. There is evidence supporting the latter hypoth- 
esis in that N III]/N V and C III]/C IV behave simi- 
larly in our overall quasar sample. In particular, when 
N III]/N V is larger/smaller than average, C III]/C IV is 
also larger/smaller by roughly the same factor. The inter- 
combination lines scale roughly together compared to the 
permitted lines of the same element. The somewhat differ- 
ent behaviors of N III]/N V and C III]/C IV in Figure 11 
(sorted by FWHM) therefore suggest that we arc under- 
estimating the weak N III] line in the broader line sources. 
The U-like shape of the N III]/0 III] and N III]/C III] 
points in this figure could then be caused by an increas- 
ing underestimate of N III] as the FWHM increases from 
^2000 km s"-'^ to ^5000 km s"-'^, offset by the increas- 
ing N/0 and N/C abundance ratios for FWHM > 5000 
km s~^. Thus, we believe this is the more likely scenario. 
This scenario may also explain the fact that the metallic- 
ities derived from all ratios agree for objects with narrow 
emission lines (see also Warner ct al. 2002: Baldwin ct al. 
2003), in which the continuum is unlikely to be systemat- 
ically overestimated. 

6. DISCUSSION 

Our main result from §5 is that the data are consis- 
tent with a trend between BLR metallicity and SMBH 
mass. The hypothesis motivating this work (§1) was that 
there should be a trend between metallicity and SMBH 
mass with roughly the same slope as the mass-mctallicity 
relation in elliptical galaxies, if the BLR gas was en- 
riched by stars in the AGN's host galaxy. We use a 
plot of Mg2 index vs. velocity dispersion, a, in Bender 
et al. (1993) to estimate the slope of the galactic mass- 
metallicity relation. Mg2 index is a line strength parame- 
ter used to derive metallicities. Bender ct al. (1993) derive 
log(Mg2) DC 0.41 log(Z/Z0) and log(Mg2) oc 0.20 log(cr). 
Using these proportionalities, along with M oi a'^, we ob- 
tain an estimate of the slope of the mass-mctallicity rela- 
tion in log space, AZ/AM ^ 0.24. We obtain a similar 
estimate {AZ/AM ~ 0.34) using a plot of [0/H] vs. ab- 
solute blue magnitude in Zaritsky et al. (1994) combined 
with a plot of log(LB) vs. log(M) in van Albada ct al. 
(1995). We also derive log(Z) oc 0.38 ± 0.07 log(MsMBH) 
from the [Z/B] oc 0.76±0.13 logfu) relation given in Trager 
et al. (2000). We employ a x minimization routine to 
obtain linear fits to the plots of metallicity vs. SMBH 
mass shown in Figure 10. These fits yield a slope of 
AZ/AM ~ 0.2 - 0.3 for all ratios involving N V. This is 
consistent with the slope of the mass-metallicity relation 
in elliptical galaxies as expected. The slope for the average 
metallicity vs. SMBH mass is ~0.2, slightly lower than the 
slopes for the N V ratios. The slope for ratios involving 
N III] is ~0, but the slope from only the three highest mass 
spectra for these ratios is 0.1 — 0.2. This may be a lower 



bound on the slope of metallicity vs. SMBH mass because 
in the highest mass spectra, the increasing N/0 and N/C 
abundance ratios are starting to dominate the underesti- 
mation of N III] caused by the systematic overestimation 
of the continuum in spectra with broad lines. 

It should be noted that the quasar-BLR metallicities 
are several times higher than the galactic values derived 
in the studies listed above. This apparent discrepancy can 
be understood in the context of normal galaxy evolution. 
Quasar BLRs sample an earlier evolutionary stage and fre- 
quently a different physical location in galaxies compared 
to the H II regions or average spectra of stars in present 
day galaxies (Hamann & Ferland 1999). However, the 
similar slope of the mass-metallicity relationship in both 
quasars and galaxies supports the premise that the gas in 
the BLR was processed/enriched by the surrounding stel- 
lar population. 

The evidence presented here for a mass-metallicity rela- 
tionship among quasars does not prove that SMBH mass 
is the fundamental parameter affecting BLR metallicity. 
However, in these composite spectra, sorted by SMBH 
mass, the slope of the trend between luminosity and metal- 
licity is at least as steep as (and possibly marginally 
steeper than) the slope is for composite spectra sorted by 
luminosity (see Dietrich et al. 2002; Dietrich et al. 2003, 
in prep). Also, composite spectra created from different 
ranges in Eddington ratio, Lboi/Ledd, show a slight trend 
with luminosity, but no trend with either SMBH mass or 
metallicity (see Warner et al. 2003, in prep). These results 
suggest that SMBH mass could be the fundamental pa- 
rameter. Further support for mass as the fundamental pa- 
rameter comes from the similarity of the mass-mctallicity 
relationship among quasars to the corresponding relation- 
ship among galaxies. These ideas are also consistent with 
the growing evidence (e.g. Shields et al. 2002) for a close 
relationship between quasars and their host galaxies. In 
particular, the relationship between SMBH mass and host 
galaxy mass (see refs. in §1) indicates that i) SMBHs are a 
natural byproduct of galaxy formation, and ii) essentially 
all massive galaxies were at one time "active." Quasars 
at high redshifts are signposts of the last stages in the 
formation of these massive SMBHs in the cores of young 
galactic spheroids. The high metallicities near these ob- 
jects reported here and elsewhere (Constantin et al. 2001; 
Dietrich et al. 1999, 2003; Hamann 1997; Hamann & 
Ferland 1999; Hamann et al. 2002; Osmer ct al. 1994; 
Petitjean et al. 1994; Warner et al. 2002) imply that 
the gas has already (at the quasar epoch) been substan- 
tially processed/enriched by multiple generations of mas- 
sive stars. This rapid evolution at high redshifts, though, 
is well within the parameters derived in some recent sim- 
ulations, which show that the densest protogalactic con- 
densations can form stars and reach solar or higher metal- 
licities at > 6 (Gnedin & Ostriker 1997; Haiman & Loeb 
2001). 

7. SUMMARY & CONCLUSIONS 

We have investigated a large sample of 578 AGNs for 
a trend between metallicity and SMBH mass. The sam- 
ple covers a redshift range from < z < 5, six orders of 
magnitude in luminosity, and five orders of magnitude in 
SMBH mass. We estimate SMBH masses using the virial 
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theorem and formulae given in Kaspi et al. (2000). To 
improve the signal-to-noise ratio and average over object- 
to-object variations, we produce composite spectra repre- 
senting each decade in SMBH mass. Composite spectra al- 
low us to better measure weak emission lines and minimize 
the influence of irregularities in any individual quasar. Af- 
ter a power law continuum fit, multi-component Gaussian 
profiles are fit to emission lines to measure their fluxes. 
Metallicities are then estimated by comparing emission 
line ratios involving nitrogen to theoretical predictions in 
Hamann et al. (2002). Our main results are as follows. 

1) We estimate SMBH masses based on both C IV and 
Hp and find approximately a 1:1 correlation for averages 
of many objects (Fig. 5). We conclude that it is valid 
to use C IV instead of H(3 to estimate average black hole 
masses in samples of AGNs. 

2) There is no trend between SMBH mass and redshift, 
but as expected, SMBH mass shows a positive correlation 
with both luminosity and the FWHM of C IV. 

3) We observe the usual Baldwin Effect in composite 
spectra representing different ranges in SMBH mass. In 
agreement with other studies (see Dietrich et al. 2002), 
N V does not show the Baldwin Effect even though it 
seems to be present to some extent in N III]. However, 
Al III curiously does not exhibit the Baldwin Effect ei- 
ther, and if anything, its REW actually increases with 
increasing SMBH mass. 

4) We find a trend in nietallicity with SMBH mass for 
emission line ratios involving N V, but not for those in- 
volving N III] or N IV] . This is consistent with the findings 
of Dietrich et al. (2002). They find a lack of a Baldwin 
Effect for N V, which in the present study leads to emis- 
sion line ratios involving N V increasing with increasing 



luminosity. However, Dietrich et al. (2003a, 2003b) find 
that metallicities derived from ratios involving both N III] 
and N V yield results of solar or greater, for samples of 70 
high redshift (3.9 < ^ < 5.0) quasars. 

5) We conclude that the data arc consistent with a trend 
between SMBH mass and metallicity and some of the data 
(N V ratios) indicate that there is a very strong trend. 
N III] seems to indicate no trend; however, upon further 
examination, we believe that the uncertainties in N III] 
and N IV] are too large to confirm or test any trend. We 
note that the correlation between SMBH mass and FWHM 
may lead to systematic overestimates of the continuum and 
thus underestimates of N III] and N IV] in the higher mass 
spectra. 

6) We estimate the slope of the galactic mass-metallicity 
relationship in log space to be AZ/AM ~ 0.2 — 0.3 (see 
§6). This slope is consistent with linear fits to the plots 
of metallicity vs. SMBH mass in Fig. 10 that arc derived 
from ratios involving N V, as well as the plot of "aver- 
age" metallicity vs. SMBH mass. The similarity in the 
slope of the mass-metallicity relationship in both quasars 
and galaxies supports the premise that the gas in the BLR 
was processed/enriched by the surrounding stellar environ- 
ment and contributes to the growing evidence for a close 
relationship between quasars and their host galaxies. 
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APPENDIX 

We fit C IV A1549 with three Gaussian profiles. Three Gaussians are necessary to obtain a good fit because the emission 
lines are not Gaussian by nature and contain broad wings. We ignore the doublet splitting in C IV because it is small 
compared to the observed line widths. N IV] A1486, where it is measureable, is fit with a C IV profile. That is, we fit it 
with three Gaussian profiles that have their FWHMs tied to the values fit to C IV and their relative fluxes and relative 
wavelength shifts are tied to the same ratios as in the flt to C IV. For the two highest mass composite spectra, we flt only 
an upper-bound to N IV] because it is too blended with the wing of C IV to measure. He II A1640 is fit with two Gaussian 
profiles, one broad and one narrow and O III] A1665 is fit with a C IV profile. The weak and badly blended doublet Si II 
AA1527, 1533 is fit with one Gaussian profile for each component. 

Lya A1216, like C IV, is fit with three Gaussian profiles, but only the red side of the line is used for minimization 
because of contamination by the Lya forest on the blue side. Each component of the N V A1240 doublet is fit with a C IV 
profile because N V is strongly blended with the wing of Lya. This connection between C IV and N V can be justified 
because they are both high ionization lines with similar excitation properties (Ferland et al. 1998). The relative fluxes of 
the two components of the N V doublet are set to the ratio of the statistical weight, g, times the A value for each transition 
(approximately 2:1 in this case), because this yields a better fit than using a 1:1 ratio. We use atomic data obtained from 
the National Institute of Standards and Technology (http://aeldata.phy.nist.gov/PhysRefData/contents-atomic.html) to 
calculate g times the A value for each transition. The weak and blended multiplet Si II AA1260, 1264, 1265 is flt with one 
Gaussian profile for each component. 

We flt O I A1303 with a single Gaussian proflle and C II A1335 with two Gaussians, one for each component of the 
doublet. The fluxes of the two C II components arc tied to a 1:1 ratio because this yields a better fit than a 2:1 ratio. 
The FWHMs and relative wavelengths are also tied together. 

We fit the doublet Si IV AA1393, 1403 with one Gaussian profile for each of the two components. The fluxes of these 
components are tied to a 1:1 ratio because this yields a better fit than a 2:1 ratio. O IV] A1403, a blended multiplet of 
five lines, is fit with one Gaussian profile for each of its five components. The flux ratios of the components in this case 
are tied to the ratio of g times the A value for each transition using atomic data obtained from Nussbaumer & Storey 
(1982). Because Si IV and O IV] are badly blended, we tie the FWHMs and relative wavelengths together for all seven 
components of the Si IV - O IV] blend. We list the sum of the blend in Table 2. 
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We fit the N III] A1750 multiplet with five C IV profiles, one for each component of the multiplet. We use C IV profiles 
to constrain the fit because N III] is a weak intercombination line. The relative fluxes of the five components are tied to 
the ratio of g times the A value for each transition because this yields the best fit. The FWHMs and relative wavelengths 
are all tied together. 

We use two Gaussian profiles to fit C III] A1909 because two Gaussians provide a fit that is just as good as the fit 
obtained by using three Gaussians. The Al III AA1855, 1863 doublet is fit with one Gaussian profile for each component. 
The flux ratio is tied to the ratio of g times the A value for each transition and the FWHMs and relative wavelengths 
of both components are tied together. Si III] A1892 is fit with one Gaussian profile and has its FWHM and relative 
wavelength tied to Al III because it is badly blended with C III] . 

We fit a local continuum around O VI A1035, which was constrained by the flux in wavelength intervals centered at 960 
A and 1100 A. We then fit each component of the O VI A1035 doublet with two Gaussian profiles. The fiuxes are tied to a 
3:2 ratio, compromising between 1:1 and 2:1, because this yields the best fit. Where they are measureable, C III A977 and 
N III A99G are fit with one Gaussian profile each, and their FWHMs and relative wavelengths are tied together. When 
calculating fluxes from measured REWs, we scale O VI as if it were sitting on top of a segmented powerlaw continuum, 
defined by our previous continuum fit redward of 1250 A and by a powerlaw with index a = —1.76 (Telfer et al. 2002) 
blueward of 1250 A. 

Each of the five spectra is fit in this fashion, with the following exceptions for the broadest (M ~ 10^" Mq) spectrum, 
which has many extremely blended lines. In this spectrum, we tie the relative wavelengths of the three Gaussian profiles 
used to fit C IV and Lya and the two Gaussian profiles used to fit C III] to the corresponding ratios from our fit to the 
M 10^ Mq sample. We do not attempt to fit N IV] or Si II AA1527, 1533, as they are undetectable due to the broad 
wing of G IV. It should also be noted that in our fit to this spectrum, the broader Gaussian profile of our fit to He II is 
blueshifted by 15 A. 
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Fig. 1. — Normalized composite spectra are shown for all five SMBH mass bins. The spectra were normalized by a power law continuum 
fit. The horizontal dashed lines indicate the continuum level for the individual normalized composite spectra, which were vertically shifted 
for better display. The normalized continuum strength is shown for the spectrum at the bottom of the figure. This same vertical scale applies 
to all other spectra, although the tick marks for '0' and '1' are not labelled. The mean SMBH mass of each spectra is listed to the right. The 
Baldwin Effect is clearly seen in these spectra. 
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Fig. 2. — Same as Fig. 1, but with an expanded vertical scale to display the dependence of the relative strength of weaker emission lines 
as a function of SMBH mass. Strong emission lines with flat tops are truncated for easier display. 
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Fig. 3. — Multi-component Gaussian fits (smooth solid curves) to observed emission lines (jagged solid curves) for a composite spectrum 
of M ~ 10* Mq. The continuum is normalized by a power law continuum fit. The continuum (at unity) and composite line fits are plotted. 
The individual Gaussian components of the fits arc shown for the fit to N III] as dotted lines. In the other panels, the dashed and dotted 
lines represent the sum contribution of each labelled emission line. There is an unidentified bump at ~1070 A that is not attributed to O VI 
(see Hamann et al. 1998). 
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Fig. 4. — Comparison of the measured FWHM of C IV to that of if/3 for all objects in our sample with observations of both lines. In five 
cases, obvious narrow Hf3 emission was subtracted prior to measuring its FWHM. The open boxes represent five composite spectra created 
from objects in selected rajiges of SMBH mass. The error bars drawn on the boxes represent the ± Icr ranges for the FWHMs of C IV and if/3 
for these composite spectra. The filled diamonds represent averages of objects in selected ranges of C IV FWHM (< 2000 km/s, 2000-3000 
km/s, 3000-4000 km/s, and 4000-6000 km/s). The dotted line is a 1:1 ratio and the dashed line is the expected ratio of FWHM(C IV) = V2 
FWHM(/i/3). 



13 



> 
O 

o 



CO 
O 



m 

cn 
O 




5 



6 



7 8 9 

_og BH Mass from Hj3 



10 



1 1 



Fig. 5. — Comparison of the SMBH mass derived from C IV to that derived from Hfd for all objects in our sample with observations of 
both lines. The filled diamonds represent averages of objects in selected SMBH mass ranges. The dotted line is a 1:1 ratio (the RMS scatter 
about the 1:1 relationship is 0.39 for the individual objects and 0.17 for the averages). For composite spectra of many objects, both methods 
yield approximately the same mass. 
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Fig. 6. — Redshift distribution of the AGN sample as a function of SMBH mass. The filled diamonds represent radio-quiet quasars and the 
open diamonds represent radio- loud quasars. There is clearly no trend in SMBH mass with redshift (for z > 0.3). 
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Fig. 7. — Continuum luminosity versus SMBH mass. There is a very strong correlation with little scatter. Symbols as in Fig. 6. 
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Fig. 8. — The FWHM of C IV versus SMBH mass. A trend is seen, but with a large amount of scatter. Symbols as in Fig. 6. 
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Fig. 9. Emission line equivalent widths as a function of SMBH mass. We calculated linear fits to the EWs for the entire SMBH mass 
range using a chi-square minimization routine (dotted lines). The slopes of these fits are given in the upper right hand corner of each plot. 
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Fig. 10. — Metallicities derived from comparisons of different line ratio diagnostics to theoretical results from Figure 5 in Hamann et al. 
(2002) are shown as a function of SMBH mass. The average includes N III]/C III], N III]/0 III], N V/C IV, and N V/O VI. 



19 



0.5 



0.4 



^ 0.3 
cr 

c 0.2 



0.1 

0.0 
1 .2 

1 .0 

0.8 

0.6 

0.4 



C III] / C IV 



D 



c 



<> 



<> 



N III] / III] 



<> 



<> 



<> 



<> 



<> 



0.6 " 

0.4 - 

0.2 - 

0.0 _ 

1 .5 " 

1 .0 - 

0.5 - 
0.0 _ 



N III] / N V 



<> 



N V / C IV 



<> 



<> 



N III] / C III] 



0.6 
0.4 

0.2 

0.0 
1 2 

10 
8 
6 
4 



<> 



N V / He 



2 - 
_ 



<> 



<> 



<> 



<> 



2000 4000 6000 8000 2000 4000 6000 8000 2000 4000 6000 8000 
FWHM FWHM FWHM 



Fig. 11. — Emission line ratios are shown as a function of the FWHM of C IV. 



Table 2 



Emission Line Data 



Line 




- 1O''M0 - 






- IO^'Mq - 






- 1O*M0 - 
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Flux/LycK 


REW" 


fwhm'' 
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0.023 
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27 


5800 
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24 


4900 
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19 
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16 


8600 


Lya A1216 
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42 
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4 
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2 
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2 
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16 
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12 
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11 
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10 
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2" 
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3100 


0.597 
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41 
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28 


6500 


He II A1640 
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13 
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13 


4200 
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7 
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3 
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8 
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9 
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5 
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4 
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4 
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3 
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^In units of A 

^^In units of km s""*" 
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